Many-Body Force and Mobility Measurements in Colloidal Systems 
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We demonstrate a technique for simultaneously measuring each component of the force vectors 
and mobility tensor of a small collection of colloidal particles based on observing a set of particle 
trajectories. For a few-body system of micron-sized polymer beads in oil separated by several 
particle radii, we find that the mobility tensor is well-described by a pairwise Stokeslet model. This 
stands in contrast to the electrostatic interactions, which were found to deviate significantly from a 
pairwise model. The measurement technique presented here should be simple to extend to systems 
of heterogeneous, non-spherical particles arranged in arbitrary 3D geometries. 



Typical approaches to understanding the macroscopic 
behavior of colloidal materials are based on microscopic 
models including thermal fluctuations and the interac- 
tions of pairs of particles. However, collective behav- 
iors of many particles are fundamentally multi-coordinate 
phenomena that can only be reduced to pair interactions 
in special circumstances. To that end, techniques for 
measuring forces in multi-coordinate systems are needed. 

Most force measurements are based on Hooke's law; 
forces of interest are inferred by measuring the deflec- 
tion of carefully calibrated springs. In the microscopic 
domain, the surface forces apparatus (SFA) and atomic 
force microscope (AFM) use mechanical cantilevers [lj. 
Static optical tweezers (OT) measurements use an op- 
tical potential. These techniques allow sensitive mea- 
surements of interactions between pairs of surfaces, but 
SFA and AFM are all but impossible to extend to multi- 
particle force measurements. Static optical tweezers can 
be used to simultaneously measure forces between sets 
of particles @, 0] , but such measurements require care- 
ful calibration of many traps and may be confounded 
by optical interactions between particles induced by the 
tweezers themselves. 

Hooke's law is not the only way to measure a force; 
forces may also be inferred from the trajectories of mov- 
ing particles. If we wanted to measure the force of grav- 
ity on a baseball, we could weigh it on a scale, but 
we could also toss it in the air and measure its trajec- 
tory. Astronomers and particle physicists infer forces 
from trajectories because appropriate calibrated springs 
are not available. Such dynamical force measurements 
have the advantage of not requiring any carefully cal- 
ibrated springs and avoid perturbing the system being 
studied. Crocker and Grier introduced a dynamical inter- 
action measurement suitable for colloidal systems called 
Markovian Dynamics Extrapolation (MDE) In this 
technique, data are collected by first positioning parti- 
cles in some initial configuration using optical tweezers, 
and then switching the tweezers off and recording the 



trajectories of the freely interacting particles using video 
microscopy. 

In Ref. 0, we introduced an alternative technique for 
extracting forces from similar trajectory data. Whereas 
the original MDE analysis requires the entire configu- 
ration space to be sampled before the interactions in 
any particular configuration can be inferred, our anal- 
ysis uses only local information to calculate forces. This 
advantage becomes crucial when measuring interactions 
in multi-particle systems because the size of the con- 
figuration space explodes as the number of particles is 
increased. We have recently used this technique to in- 
vestigate surfactant-controlled electrostatic interactions 
between pairs of colloidal particles in a non-polar sol- 
vent [f|, and deviations from pairwise additivity for the 
forces between sets of several particles in the same system 
Q • These measurements have implicitly involved several 
coordinates — even a system of two spheres moving in two 
dimensions has a total of four coordinates — but for ease 
of presentation, we have previously projected the mea- 
sured forces onto a single coordinate. Here, we describe 
the extension of the technique to multiple coordinates 
and exhibit the full force and mobility curves for each 
component of a multi-particle system. 

The central result of Ref. [5( is a relation between the 
force, /, and the two statistical parameters that describe 
ensembles of short-time Brownian trajectories, the drift 
velocity, Vd, and the diffusion coefficient, D: 



f = k B TD- l v d . 



(1) 



This is the correct result for a single particle in ID sub- 
ject to an external field. As demonstrated in Ref. [f|, 
this scalar equation also holds separately for each of the 
hydrodynamic normal modes of a multi-coordinate sys- 
tem, with each normal mode having it's own diffusion 
coefficient. 

To derive this relation, we start with the expression for 
the drift velocity of a particle in a heavily over-damped 



viscous medium: 



v<i = bf. 



(2) 



where b is the particle mobility. If the particle mobility 
were already known, then this relation would be suffi- 
cient to determine the forces from particle trajectories 
simply by measuring Vd- In principle, the mobility can 
be calculated from the properties of the particles and the 
solvent. Alternatively, the mobility can be eliminated 
from Eq. (JSJ) using the Einstein relation: 



b = (3D. 



(3) 



A simple rearrangement then gives Eq. fl} . Since D can 
be measured directly from the particle trajectories, this 
eliminates the need to know the mobility ahead of time. 

This approach is readily generalized to a system of 
many particles in multiple dimensions. In that case, hy- 
drodynamic interactions lead the drift velocity of one par- 
ticle to be coupled to the force on all the others. The drift 
velocity and force vectors, and f , are linearly related 
through the mobility tensor, b: 



v d = bf . 



(4) 



The vector drift velocity and force are made up 
of the concatenation of the components of the ve- 
locity and force for each particle, e.g. = 
{Vd y n } an( i similarly for 

the force vector. Hydrodynamic interactions lead to off- 
diagonal terms in the mobility tensor, so that a force on 
one particle may produce a velocity on another particle, 
or a force directed along one coordinate of a particle may 
produce a velocity along a different coordinate Q . 

A multi-coordinate Einstein relation holds for the mo- 
bility and diffusion tensors 0, [13] : 



b = /3D, 



(5) 



so we can write down a multi-coordinate version of Eq. 
© : 



f - k B TT)-\ d . 



(6) 



Since the mobility tensor is symmetric, it can be diag- 
onalized by a set of normal modes. The scalar Eq. ((T|) 
holds separately for the force on each normal mode of the 
mobility tensor, as was shown in Ref. [j|. However, it 
isn't necessary to solve for the normal modes in order to 
calculate each component of the force on a multi-particle 
system. Eq. © can be evaluated directly. Contrary to 
the assumptions of Ref. [Bj] , there is no need for the force 
to be directed along a normal mode, nor is there any 
restriction to conservative forces. 

Over appropriate time scales, the mean displacement 
and displacement covariance increase linearly in time, 
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FIG. 1: Mean Displacement and Displacement Covariance (a) 
Mean displacement and (b) displacement covariance versus 
time for each coordinate of three particles arranged in an equi- 
lateral triangle as shown with side length s — 4.4a = 2.6 //m. 
Red lines are best fits of Eqns. (FT} and ||8} to the data. 



with slopes related to the components of the drift ve- 
locity, Vd,i, and diffusion tensor, Dij [ill ]: 

(xi(t + r) - Xi(T)) T = (v dti + ^2 d kD k i)t, (7) 



CaV T (Xi(t+T)-Xi(T),Xj(t+r)-Xj(T)) = ■ r,,. (8) 

A subtle point here is that the increase of the mean dis- 
placement with time given by Eq. ([7|) is due not only 
to the drift velocity, but also depends on the local gradi- 
ent of the diffusion tensor. This correction will become 
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FIG. 2: Velocity and Diffusion (a) Drift velocity vector and 
(b) diffusion/mobility tensor for particles arranged in an equi- 
lateral triangle as a function of side length, s. Diffusion 
(mobility) values are normalized to Do = ksT /6nr]a = 117 
nm 2 ms _1 (6o = l/Qivqa = 29.5 j ums _1 pN _1 ). Red lines on 
the diffusion/mobility tensor plot are predictions based on 
Eqns. Ull) and (fT2)l . 



important when particle surfaces are sufficiently close to 
each other, but it is negligible in the data presented here. 

This suggests a straightforward approach to measur- 
ing the force and mobility in multi-particle systems. We 
measure the velocity and diffusion by fitting these lines 
to the statistics of observed trajectories, taking averages 
for a given delay time, t, over the set of non-overlapping 
pairs of sample times separated by t. The intercept term 
in the displacement covariance equation, , arises due to 



a combination of uncertainty in position measurements 
and finite shutter speed 

In the following, we apply this technique to measure 
the force and mobility for sets of three particles arranged 
in an equilateral triangle and for seven particles arranged 
in a hexagon with a particle at its center. The particles 
are poly(methyl methacrylate) spheres of radius a = 600 
nm suspended in a solution of 500 (aM NaAOT in hex- 
adecane 

The mean displacement and displacement covariance 
versus time corresponding to Eqns. ([7]) and |(8]) are shown 
in Fig. [1] for particles arranged in an equilateral triangle 
with side length s = AAa = 2.6 jirsx. This procedure 
is repeated at a series of side lengths, s. The measured 
drift velocity and diffusion tensor are shown in Fig. [2] for 
particles arranged in an equilateral triangle with a series 
of side lengths, s. 

With both Vd and D in hand, calculating the entire 
force vector is a simple application of Eq. ([5]). The full 
force vector is shown for a series of separations in Fig. [3] 
for both the equilateral and hexagonal particle configu- 
rations. These are the same data as in Ref. Q, where 
we reported only forces on the breathing mode calculated 
using Eq. |T]). 

There is nothing that limits this analysis to 2D. If 3D 
position data were available, it could be taken into ac- 
count with no additional complications. In fact, it would 
generally be necessary to use 3D data in order to correctly 
infer the interparticle forces (l3| : however, if all the par- 
ticles are in plane and the surroundings are isotropic, the 
force measurement splits as 



k B TT) x ^v dtXy , 



(9) 
(10) 



We expect the forces to be central and thus in plane, so 
we need only observe the in plane position in order to 
find the forces. The observed trajectories must be short 
enough that the particles remain roughly in plane. We 
ensure that the surroundings are effectively isotropic by 
conducting our experiments more than 30 particle diam- 
eters from any wall. 

To this point, we have refrained from giving any inter- 
pretation for the form of the velocity, diffusion, and force 
versus separation curves in order to emphasize that our 
force calculation is independent of any model of the par- 
ticle interactions. Our only assumptions are that 1) the 
Einstein relation, Eq. ((5]), holds; 2) the drift velocity and 
diffusion coefficients may be considered locally constant 
over a single trajectory; and 3) the gradient of the diffu- 
sion/mobility tensor can be neglected in determining the 
drift velocity. We satisfy conditions 1 and 2 by choos- 
ing appropriate experimental timescales, i.e. the frame 
rate and number of frames in single trajectory from Eq. 
(J?]). We can check condition 3 by examining the diffu- 
sion tensor data in Fig. [51 The magnitude of V ■ D is 
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FIG. 3: Forces versus Separation Components of the force as a function of side length for (a) each of three particles arranged 
in an equilateral triangle with side length s and (b) seven particles arranged in a hexagon with a particle at its center. 



less than 0.1 /im/s at the smallest measured separation, 
which is a negligible fraction of the drift velocity. No- 
tably, we have not made any assumption that the forces 
or hydrodynamic coupling are pairwise. With the present 
technique, we can measure whether the interactions are 
pairwise. 

Having established the general force measurement 
technique, we now explore the structure of the multi- 
particle diffusion/mobility tensor. The interparticle 
forces in this data set were interpreted in Ref. [7|, 
where we found good agreement with linearized Poisson- 
Boltzmann theory so long as particle surfaces were taken 
to be at constant electrostatic potential rather than con- 
stant charge density. These boundary conditions lead to 
the observed non-pairwise additivity of the forces. 

In these data, where the separation between particles is 
always greater than 4 particle radii, the mobility data are 
well-described by the Stokeslet Superposition Approxi- 
mation (SSA) [141 ]. In this approximation, the mobility 
is expanded in powers of 1/r and truncated at first order. 
More accurate mobility calculations for smaller separa- 



tions are described in Ref. [151 ]. 

The mobility tensor can be split up into particle- 
particle blocks (2 x 2 in this 2D case). The diagonal 
blocks represent the self-mobility of a single particle and 
are given in the SSA by 

b S = -^—12 = W2, (11) 

bTrr/a 

where rj is the dynamic viscosity of the fluid, a is the 
radius of the particle, and I2 is the 2x2 identity. The 
off-diagonal blocks represent hydrodynamic interactions 
between a pair of particles, i.e. the cross-mobility, and 
are given in the SSA by 

where r is the separation between the particles, and ff 
is a projection operator that projects the force along the 
pair separation axis. In the fixed camera basis, its com- 
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FIG. 4: Self-Mobility and Cross-Mobility Self blocks (a,c) and cross blocks (b,d) of the diffusion/mobility tensor for equilateral 
(a,b) and hexagonal (c,d) particle configurations. Red lines are the predictions of Eqns. Qllj l and (|14p . The self-mobility blocks 
are expressed in fixed camera coordinates, and the cross-mobility blocks are expressed in pair coordinates as described in the 
text. The self-mobility blocks reference only a single particle, so the horizontal axis is labeled by the side length, s, of the 
configuration. The cross-mobility blocks reference two particles, so the horizontal axis is labeled by the separation, r, between 
the pair. In the equilateral configuration, all particles are separated by the same distance, so there is no distinction between 
r and s, but in the hexagonal configuration, there are several different pair separations, r, for a single side length, s. In the 
cross-mobility plots (b,d), each point has a line through it indicating the direction of the pair separation axis in fixed camera 
coordinates. 



ponents are 



cos 2 ((9) sin(0)cos(0) 
sin(0)cos(0) sin 2 (0) 



(13) 



where is the angle between the pair separation axis and 
the camera x-axis. Predictions of this theory based on 
the nominal particle radius, a = 600 nm, are shown as 
red lines in Fig. [21 and are in good agreement with the 
measurement. 

In order to further bring out the structure of the diffu- 
sion/mobility tensor, we note that for identical particles, 
in this approximation the self-mobility blocks should all 



be identical and constant, and the cross-mobility blocks 
should all collapse when expressed in coordinates adapted 
to each pair with basis vectors parallel and perpendicular 
to the separation axis (0 = 0). In these pair coordinates, 
the components of the cross-mobility blocks are simply 



1 



8Trrjr 



2 
1 



(14) 



We plot measured values of the self-mobility blocks as 
a function of side length, s, and the cross-mobility blocks 
as a function of separation, r, for both the equilateral 
and hexagonal configurations in Fig. 3J For the cross- 
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FIG. 5: Normal mode mobility and projected force (a) Mobility versus side length for each of the 6 linearly independent normal 
modes of the equilateral configuration, (b) Force projected onto each of the normal modes. Only the breathing mode has a 
large non-zero force. Black points with error bars show the result of analyzing the perimeter coordinate with the scalar Eq. 
©, as in Ref. Q. 



mobility blocks, the direction of the interparticle separa- 
tion is represented by the direction of a blue line through 
each point. The Stokeslet predictions are plotted as red 
lines. Agreement is generally quite good, although the 
diffusion coefficient along the camera's ir-axis is slightly 
larger than the diffusion coefficient along the y-axis. We 
believe that this is due to anisotropic fluctuations in the 
microscope stage, which lead to a small increase in the 
apparent diffusion coefficient along the x-axis. This effect 
is also visible in the cross-mobility blocks, where a small 
systematic dependence on the direction of the pair sep- 
aration axis can be discerned. The consistent collapse 
of these curves demonstrates that hydrodynamic inter- 
actions are essentially pairwise for these separations, in 
contrast to the electrostatic interactions in this system 
which showed significant departures from the pairwise 
model [2j. 

Analysis of the normal modes of the diffusion/mobility 
tensor gives additional insight. From the mobility point 
of view, when a force is applied along one of these normal 
modes, it induces a velocity in that normal mode only and 
does not couple to the other normal modes. From the dif- 
fusion point of view, fluctuations in these coordinates are 
uncorrelated. Fig. [5] shows the mobility as a function of 
side length of the equilateral triangle for each of the 6 lin- 
early independent normal modes |16j . Two of the normal 
modes are unique; the others pair to form two degenerate 
subspaces. The degeneracy is a consequence of the 3-fold 
rotation symmetry of the equilateral configuration. Ro- 
tating any normal mode by 27r/3 leads to another normal 
mode with the same eigenvalue. The original mode and 
the rotated mode will then span a degenerate subspace 



unless they are in fact the same mode — that is, unless the 
original mode possesses the full 3-fold symmetry of the 
particle configuration. Fig. [5b shows the magnitude of 
the projection of the total force onto each of the normal 
modes. Since the particles are nearly identical and the 
force is central, the force is directed almost entirely along 
the breathing mode. 

Given this observation, it is interesting to consider 
whether the force could be characterized by reducing 
the system to ID at the outset by considering only the 
perimeter of the equilateral triangle at each time and 
calculating the force using the scalar Eq. (TT]). An advan- 
tage of this analysis is that the perimeter is insensitive 
to drifts and fluctuations in the microscope stage. The 
result of this procedure is shown in Fig JS]as black points 
with error bars, as presented in Ref. 7]. In these data, 
there is good agreement with the forces calculated by 
projecting at the end, but such agreement should not 
be taken for granted. When working with non-Cartesian 
coordinates (like the perimeter), it is generally necessary 
to make geometrical corrections to Eq. (J7|) that reflect 
variations in the volume of phase space associated with 
a given coordinate. The necessary corrections can be de- 
rived by expanding the mean and mean square displace- 
ments of the perimeter (or other coordinate function) 
in small displacements around the initial configuration, 
keeping terms only up to first order in dt. 

As an example, consider the case of two particles. If 
there are no forces, then according to Eq. (j4|) there 
should be no drift velocity for any of the Cartesian coor- 
dinates of the particles. Nevertheless, we expect diffusion 
will cause the mean separation between the particles to 



7 



increase in time. If we consider the separation as a single 
coordinate and naively equate the slope of the mean dis- 
placement versus time with the drift velocity, then upon 
applying Eq. |T]) we will find a non-zero "effective force." 
For the separation, r, between two particles, the effective 
force is f c g = ksT/r. For particles separated by 2 mi- 
crons at room temperature, the effective force is 2 fN, 
which is just below our experimental resolution. 

Measurements using this technique have already pro- 
vided valuable insight into the charging mechanism of 
colloids in non-polar solvents. Soon, we hope to examine 
systems of particles with dissimilar charges and aqueous 
systems. Such force measurements are not limited to 
translational coordinates. It should be possible to use 
the same technique to measure torques on anisotropic 
particles. 
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